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ABSTRACT 

We present a simple resonant Hamiltonian model for the vertical response of a stellar disk to the 
growth of a bar perturbation. As a bar perturbation grows stars become trapped in vertical Inner 
Lindblad resonances and are lifted into higher amplitude orbits. The vertical structure of a boxy and 
peanut shaped bulge as a function of radius and azimuthal angle in the galaxy plane can be predicted 
from the strength and speed of the bar perturbation and the derivatives of the gravitational potential. 
This model predicts that stars on the outer side of the resonance are lifted higher than stars on the inner 
side, offering an explanation for the sharp outer edge of the boxy/peanut. 

Subject headings: 



1. INTRODUCTION 

The existence of spiral galaxy bulges with boxy or 
peanut shapes has been recognized for many years through 
the study of edge-on galaxies (e.g., Burbridge & Burbridge 
1959; de Vaucouleurs 1974). Recent work has found that 
edge-on galaxies with prominent boxy or peanuts shaped 
bulges in most cases show kinematic evidence for non- 
circular motions associated with a bar, whereas galax- 
ies lacking boxy or peanut shaped bulges do not show 
kinematic evidence for a bar (Bureau & Freeman 1999; 
Merrifield & Kuijken 1999). A few cases exist where 
highly inclined systems exhibit both bars and peanuts 
(e.g., NGC 7582, Quillen et al. 1997 and NGC 4442, Bet- 
toni & Galletta 1994). The link between a galactic bar and 
a boxy/peanut shaped bulge has been firmly established. 

Mechanisms proposed to explain the formation of boxy 
and peanut shaped structures include accretion of small 
satellite galaxies (Binney & Petrou 1985), bar-buckling 
via the fire-hose type instability (Toomre 1966; Raha et 
al. 1991; Merritt & Sellwood 1994; Fridman & Polyachenko 
1984) and resonant heating (Pfenniger 1984, 1985; Combes 
et al. 1990). The fire-hose (bar-buckling) instability is 
a global time dependent instability which relies on the 
centrifugal force caused by highly eccentric orbits in the 
plane caused by the bar and is damped via a dissipa- 
tive process similar to Landau damping (Fridman & Poly- 
achenko 1984). The resonant heating model has been 
inspired by the study and classification of orbit families 
in 3-dimensional systems (Pfenniger & Friedli 1991) and 
through the phenomena displayed in N-body simulations 
(Combes et al. 1990). 

No study has yet considered the evolution of 3 dimen- 
sional stellar orbits in a barred system that is varying with 
time. In this paper we investigate the possibility that the 
growth of the bar itself traps particles in the vertical Inner 
Lindblad Resonances (ILRs). We develop a simple Hamil- 
tonian analytical model for these resonances and use it 
to predict the vertical distribution of a bar following the 
growth of the bar. Our approach of considering resonance 
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trapping is similar to that explored by Sridhar & Touma 
(1996) for disk heating. However their model is for the 
general case of a growing axisymmetric disk, whereas our 
model applies specifically to a growing bar perturbation. 

2. DYNAMICAL MODEL 

To exihibit resonant phenomena Contopoulos (1975) 
showed that the Hamiltonian must be expanded to at least 
third order in the epicyclic approximation. When this is 
done there are higher order cross terms in the action vari- 
ables and the model will exhibit phenomenology such as 
the bifurcation of orbit families seen in characteristic dia- 
grams which are created through numerical integration of 
orbits. 

We first consider the Hamiltonian lacking non- 
axisymmetric or time dependent perturbations 



H (L,p r ,p z ;9,r,z) 



P 



+ ^ + ^ + Vo(r,z) (1) 
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where Vq is the gravitational potential which is assumed to 
be axisymmetric and L,p r ,p z are the momenta conjugate 
to the cylindrical coordinates 9, r, z. For the 2 dimensional 
problem restricted to the Galactic plane (z — 0,p z = 0), 
Contopoulos (1975) showed how to put the Hamiltonian 
in the form 

K (h,I 2 ;9i,h) = h + Kl l + ni 2 +all + 2bI 1 I 2 +cll.. (2) 

in a third order epicyclic approximation. The action vari- 
ables I\ = y~ J rdr and 1% = Jo — J c are integrals of 
motion. Jo is the particle's angular momentum, and J c 
is the angular momentum of a particle in a circular orbit 
at a radius r c which is the radius of a circular orbit with 
energy h. 9i is the epicyclic angle and 9 2 is the azimuth of 
the epicyclic center. 9\ and 9 2 are the angle variables con- 
jugate to I\ and I 2 . k is the epicyclic frequency and £1 the 
angular rotation rate; both are evaluated at r — r c , z = 0. 

To extend the theory into the third dimension we must 
add additional terms to the restricted Hamiltonian which 
depend upon z; 

H =K Q + H 0z . (3) 
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We expand the potential Vq about z and r c defining the 
variables 

2 d 2 V 



A = 



dz 2 
d^Vp 
<9z 4 



(4) 



(5) 



which are evaluated at z = 0, r = r c . Inserting this into 
the full Hamiltonian we find 

H 0z ( L, Pr , Pz ; 0, r, z) = \ (p 2 z + u 2 z 2 ) + 

+£f [r c fl 2 (r - r c ) + ±(« 2 - 3fl 2 )(r - r c f] 
where we have made use of the relations f2 2 = f^?, 

k 2 = + 3S1 2 and both £1 and k are evaluated at 

r = r c , z = 0. 

For this Hamiltonian following Contopoulos (1975) we 
can choose a third action variable I 3 to represent the am- 
plitude of vertical oscillations. 



cos(0 3 ) 



(6) 



^«-(2/ 3 ^) 1 / 2 sin(0 3 ) 



where v is the vertical oscillation frequency and 3 is the 
angle associated with the vertical oscillations. The above 
expressions are given to first order in l\J 2 and require 
higher order terms to be exact. 

After a couple of canonical transformations (see ap- 
pendix) the Hamiltonian to second order in I^ 2 takes the 
form 



H 0z (I u I 2 ,I 3 ; 0i, & , 3 ) = vh + 3^ + K - 3ft 2 ) ^ 



, f„2 or>2\ vI lh 

Ak 
(7) 

Since this is independent of the angles 9i,9 2 ,9 3 , the mo- 
menta Ii, I 2l I3 are conserved quantities and we say we 
have put the non-perturbed Hamiltonian, (Eqn. 5), in ac- 
tion angle form. 

2.1. Perturbation 

We now consider the form of the perturbation of the 
gravitational potential caused by a bar. In the plane of 
the galaxy we can expand the gravitational potential in 
terms of Fourier components 

(8) 



Vi(r, 0, z = 0) = A ™ cos(m(0 - Sl b t)) 



(for example see Quillen et al. 1994). The strongest term 
due to the bar should be the m = 2 term with A 2 < so 
that the bar major axis lies along — fltf = 0, n. For small 
z we assume that the potential can be expanded so that 
we can write for each m Fourier component 

V lm (r, 0, z) = A m cos(m(0 - fi 6 i)) (\ - (9) 

where h\ is the scale height of the perturbation. We also 
assume that A m and hi vary only weakly with radius. 
Since z 2 2i 3 -cos 2 (0 3 ) we can rewrite the perturbation 

components to second order in I^ 2 
V lm (h,6 3 ;I 2 ,e 2 ) = cos(m(0 2 - Q b t) ± 20 3 ) + ... 



Ah 2 v 



We focus on the particular resonance that we expect is 
important and ignore all the other terms. 

H( I 2 ,h; 2 , 3 ) - vh + + O/2 

+/ 3 ^cos(m(02-r! fc t)- 20 3 ) + ... (11) 

To go into a coordinate system with an angle that li- 
brates when a particle is trapped in the resonance we do 
a canonical transformation with the following generating 
function 

F 2 (J 3 , J 2 ;0 3 ,0 2 ) = J2O2 - y[m(0 2 - n b t) - 20 3 ] (12) 
so that our resonant angle and new momenta 
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<P = e 3 --(e 2 -n b t) 

1-2 — <H — J3~2 



(13) 



h = J 3 

and 02 is unchanged. Our generating function is time 
dependent; dF 2 /dt = ^£L b J 3 . Following this canonical 
transformation our Hamiltonian becomes 



h( j 2 , j 3 ; 2 , (f>) = [v- f (n - n 6 )] j 3 + 

+J 3I ^cos(2^) + OJ 2 + ... 



32i/ ; 



j2 



(14) 



The motion in the two different degrees of freedom are now 
decoupled so we can neglect J 2 and write the Hamiltonian 
in the simpler form 

H{J 3 ;4>) = z/ J 3 + ajj - e,J 3 cos(2</>) 



where 



f(0 



fit) 



(15) 



(16) 



A 



(10) 



32j/ 2 

A 



Equation(15) has resonant term proportional to J 3 in- 
stead of proportional J 3 ^ 2 as was the case for the inner 
Lindblad resonance, considered by Contopoulos (1975), 
and for e-e' orbit-orbit resonances commonly studied in 
Celestial mechanics (Murray & Dermott 1999). However 
the Hamiltonian is similar to that of J+2:J orbital reso- 
nances studied by Borderies & Goldrcich (1984). Orbits 
captured in the 2:1 vertical ILR were seen in the numerical 

simulations of Combes et al. (1990) however they expected 

1/2 

the resonance to be first order in J 3 ' rather than second 
order as we have found here. By symmetry the gravita- 
tional potential when expanded should not contain terms 

1 12 

proportional to z, and so terms oc J 3 , unless we consider 
a global instability such as the fire-hose bar buckling in- 
stability. Indeed Pfenniger & Friedli (1991) found that a 
temporal z-symmetry breaking (due to a collective effect) 
would make the resonance first order. Friedli & Pfenniger 
(1990) reported that N-body simulations which forced z- 
symmetry did grow boxy/peanut shaped bulges but more 
slowly than those which did not force the symmetry 

Unless the galaxy is extremely thin and the bulge small 
we expect that the 2:1 vertical ILR will occur in the in- 
ner region of the bar (Pfenniger & Friedli 1991; Combes 
et al. 1990). It is also possible that the 1:1 vertical ILR 
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will be present. The 2:1 vertical ILR produces banana- 
shaped orbits, referred to as BAN in the characteristic di- 
agrams of Pfcnniger & Friedli (1991) and the 1:1 vertical 
ILR produces titled orbits that are referred to as anoma- 
lous (ANO) by Pfenniger & Friedli (1991). If the galaxy 
is very thin it is also possible that there will be orbits 
associated with the 4:1 vertical ILR (Pfenniger 1985). 

If we concentrate on the 2:1 vertical ILR then the per- 
turbation term from the bar that is relevant is the m = 4 
Fourier term; if we focus on the 1:1 vertical ILR then the 
m = 2 term is relevant. It has long been known that bars 
exhibit significant m = 4 Fourier components and contain 
boxy or square contours when projected onto the plane of 
the galaxy (Athanassoula 1990). From Figure 4 by Quillen 
et al. (1994) we estimate that the m = 4 Fourier compo- 
nent of the barred galaxy NGC 4314 the potential peaks at 
A4 <~ —0.04 in units of v 2 and at A 2 <~ —0.16 in the same 
units. We can use these as typical values for the potential 
perturbations in a barred galaxy. 

Fig. 1- Contour levels of the Hamiltonian in 
Equation(17) have been plotted where x = \/2j3 cos cfi 
and y = sin <j>. Only one resonance region exists 

for ^ > 1, two exist for — 1 < ^ < 1 and three for 
^ < — 1. The resonance bifurcates at ^ = ±1 giv- 
ing banana-shaped orbits at <p = and upside down 
banana-shaped orbits at <f> = it. The bar strength 
sets e and v varies with radius from the center of the 
resonance where v = 0. Between '4 = — 1 — > 1 no 
small js circulating orbits exist. We can think of the 
progression of resonance sizes between " ~ — 1 — > 1 
as describing the radial variation of the width of the 
peanut. We expect maximum vertical amplitudes at 
— ~ — 1 which is on the outer, larger radius side of the 
resonance. As the bar grows phase space moves down- 
wards from plot to plot and particles are trapped and 
then lifted upwards. Likewise if the bar slows down 
the same thing will happen. 



at 14 ± 4, and the Milky Way bar likely to be lower than 
the typical value. Using these typical values we estimate 
that h\/r res <~ 0.3. 



-3- 



0.08 



h\ I v res 
0.3 

■An/Vc 

0.04 



(19) 



We expect the resonance to be important in the region 
where v — — > e or across a range dr; 



dr 



~A r 



v Av 2 h\ 



(20) 



We now transform Equation (15) so that our coefficients 
are unitless. We put length in units of r reSl the cen- 
ter of the resonance where v' = and time in units of 
l/v. Our action variables are in units of r 2 es v so we de- 
fine j'3 = Js/r 2 es v. We divide the entire Hamiltonian by 
r 2 es v 2 so that our new Hamiltonian becomes 

Kh,4>)/{ar 2 res ) = jl + vjs - m cos(20) (17) 

where 

32z/z/ -32i/ hi 

^=^- ^ 18 

r A v r 

res res 

SA m 8A m v 2 
e = ~ — 

r 2 es hf\ v 2 v 2 r 2 es 

and we have assumed that the vertical form of the poten- 
tial can be approximated by A <~ -rz-. A should be negative 

because the vertical oscillations are expected to be slower 
for particles traveling well above the galactic plane. 

For a bar with a peak ^ <~ —0.04, as estimated from 

NGC 4314. vr res ~ 2Q,r res ~ 2v c at the the 2:1 resonance, 
and vr res <~ v c at the 1:1 resonance. We estimate that e 
peaks at ~ 0.08 for both resonances because A 2 <~ 4A4. 
Barred galaxies have an observed ratio of bar length to 
peanut length of 2.7 ±0.3 (Lutticke et al. 2000). Measure- 
ments of the ratio of bar length to thickness give a typical 
value of ~ 10 (Sellwood & Wilkinson 1993) with the mea- 
surement of Lutticke et al. (2000) being on the high end 



where jr denotes a derivative with respect to radius. If v 
and v 1 varies slowly with radius then the resonance can be 
felt over a significant range of radius. 




t V/S = -1 



Fig. 2- A cartoon of the edge-on galaxy labeling the 
approximate location of the peanut in terms of the 
values of ^. 
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Fig. 3- Particle distribution following growth of a per- 
turbation with e = 0.1 for different values of P. The 
initial dispersion in js was a = 0.1. The perturbation 
was grown in P = 3 periods and then the integration 
continued for 5 times longer than the growth period so 
that orbits become more evenly distributed in phase. 
x and y axes are defined as in Fig. 1. 

Given the definition of our variables we estimate that 
the expectation value 



<^ 2 } ~ (h)K 



(21) 



From Figure 1 we infer that when = = — 1 the minimum 
value of vertical oscillations occurs for the fixed points at 
J3 = 1. We therefore expect the maximum height of the 
the peanut to be at | = — 1 where < j'3 >^ e. For e <~ 0.08 
this implies a peanut width of order 2\/Zr res ~ O.Qr res 
where we gain the factor of 2 when we consider orbits of 
both signs (0 = and <fi = n). This is sufficiently large 



to explain the observed vertical extent of boxy/peanut 
shaped bulges. 

2.2. Structure of the resonance 

As shown by Borderies & Goldreich (1984) the reso- 
nance contains 3 separate regions. To show the separate 
regions we use a coordinate system x = y/2j3 cos <ft and 
y = V2j3 sin 0. Contours of constant h (Eqn. 17) are 
shown in Figure 1. For j < — 1 three fixed points exist, 
for —1 < I < 1 two fixed points exist and for j > 1 only 
one fixed point exists. Orbits circulating around = 
correspond to banana-shaped orbits with positive z at the 
ends of the bar, and orbits correspond to banana-shaped 
orbits with negative z at the ends of the bar. As shown 
in Figure 8b of Pfenniger & Friedli (1991) banana shaped 
orbits appear on the characteristic diagram at a particular 
value of the Hamiltonian, consistent with them appearing 
at a particular value of |. For orbits circulating around 
the origin, the resonant angle, 0, circulates; referred to 
as oscillation as opposed to libration in analogy with a 
pendulum. 

We consider the signs of each term in Eqns.(15,17). 
v > but we expect that A, A^ 1 A 2 < 0. Consequently 
a < and e, e > and v has the opposite sign of v' . Be- 
cause v and fl increase at small radii we expect v' > and 
v < for r > r res . We see from Fig.l that near D — —e no 
orbits exist with low values of j 3 . This implies that all or- 
bits will rise out of the plane of the galaxy. As v decreases 
the mininum size of the vertical oscillations increases until 
| < — 1, then orbits with small amplitude vertical oscil- 
lations appear again. The vertical extent of the galaxy 
should increase until a particular radius and then drop 
abruptly (sec Figure 2). This would naturally correspond 
to a bowtie or peanut shape when viewed edge-on. 

2.3. Growing the bar 

To explore the simple model given in Eqn. (17) we inte- 
grate particle trajectories assuming that the perturbation 
is growing. In other words we assume that e begins at 
zero and reaches a maximum value a few rotation periods 
later. We begin the integrations with a particle distribu- 
tion p{jz) oc exp(— j|/(2cr 2 )) and with angle 4> randomly 
chosen between and 2tt. 4000 particles were integrated 
using Eqn. (17) for the simulations shown in Figures 3-6. 
The growth time for the perturbation is given in periods 



where P = ^-\ 



Since 



-es ~ 7^32 ~ 1/3, P is 
about 1/3 the bar growth time in units of the vertical os- 
cillation period. 

In the adiabatic limit the capture probability of a re- 
gion of phase space depends upon the ratio between the 
rate that phase space volume is swept across the separa- 
trix bounding this region divided by the rate that phase 
space volume is swept across the entire separatrix bound- 
ing the entire resonance (Henrard 1982; Borderies & Gol- 
dreich 1984). For D < phase space initially looks like 
the panel on the bottom right of Fig. 1 and then moves 
upward on this figure. When the separatrix touches the 
particle distribution centered at the origin, particles are 
captured into librating regions near = and ir (see Fig- 
ures 3-4). For v > phase space initially looks like the 
panel on the top left of Figure 1 and then moves downward 
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on this figure. When v ~ e particles at low j'3 are captured 
into the librating regions around <f> = 0,7r. Figures 3-4 il- 
lustrate that smoother particle distributions are achieved 
when the growth rate is faster (shorter number of periods 
P) as we would expect when the variation is no longer 
adiabatic, and when the initial momentum distribution is 
wider (larger a). 

Given a distribution in j 3 and <fi resulting from the 
growth of the perturbation we can predict the vertical dis- 
tribution as a function of azimuthal angle in the plane of 
the galaxy. For a value of the resonant angle <f> and an 
azimuthal angle 8 — f^i in the plane of the galaxy in the 
frame rotating with the bar, #3 is determined via Eqn.(13) 
and z is then determined from the definition of our ac- 
tion angle variables (Eqn. 6) . In Figure 5 we show vertical 
distributions for the action angle distributions shown in 
Figures 3 and 4b for m = 2 and in Figure 6 we show the 
equivalent but for m = 4. 

We see from Figures 5 and 6 that the vertical height 
reaches a maximum near j ~ — 1 as we expected from the 
form of the resonance. The height increases slowly between 
I ~ 1 to — 1 and then drops swiftly, confirming concep- 
tually our bowtie picture shown in Fig. 2. For m = 2 
(shown in Figure 5) corresponding to the 1:1 vertical ILR, 
the major axis of the bar becomes bulbous, as would be 
expected from a sum of anomalous orbits of both orienta- 
tions (shown in Figure 7). For m = 4 (shown in Figure 
5) corresponding to the 2:1 vertical ILR, both the major 
and minor axes of the bar have increased vertical heights. 
Again this is expected from a sum of banana-shaped orbits 
of both orientations (see Figure 7). 

Observations of peanut shaped galaxies where both the 
bar and the peanut are observed (e.g., Quillen et al. 1997) 
suggest that the ends of the peanut are aligned with the 
bar. So it could be that the 1:1 resonance orbit family is 
more likely to explain the bulbous ends of boxy/peanut 
shaped bulges. However when the radial oscillations of 
the galaxy are taken into account the banana-shaped or- 
bits may still be viable, particularly if the extent of their 
vertical oscillation is smaller along the minor axis than it 
is along the major axis. 




Fig. 7— Anomalous orbits at the 1:1 resonance are 
shown on the top and give increased widths at the 
ends of the bar. Banana shaped orbits at the 2:1 res- 
onance are shown on the bottom and give increased 
widths at the ends and sides of the bar. The x axis is 
assumed to lie along the major axis of the bar. 

We remind the reader that in Section 2.1 we dropped 
the last term in Eqn. (7) which is proportional to As 
a bar grows the epicyclic amplitude I\ of the stars will 



grow. This has the effect of shifting the location of the 
vertical resonance so that v 1 becomes 

v ' = v -—(£l-n b ) + U 2 -SO 1 ) — . (22) 
2 4k 

For a flat rotation curve n 2 = 2il 2 so we expect that v' 
will be more negative than expected if this additional term 
is neglected. This has the affect of moving the resonance 
outward in radius. One reason that the 2:1 ILR resonance 
was considered a logical choice for exciting the peanut is 
that when the galaxy is thin the 1:1 may only occur at 
extremely small radii, well inside the galaxy bulge. How- 
ever the estimates for the resonance locations were made 
assuming 1\ ~ (e.g., Pfcnniger & Norman 1990). The 
point here is that the growth of the epicyclic amplitude 
will shift the location of the vertical resonances outwards 
so that the lower order resonances may manifest at larger 
radii than previously considered. It is possible that the 
growth of the bar causes the vertical resonances to sweep 
through the inner part of the galaxy. 

3. SUMMARY AND DISCUSSION 

We have presented a simple Hamiltonian dynamical 
model for the vertical inner Lindblad resonances in the 
presence of a bar perturbation. The Hamiltonian model 
is similar to that introduced by Contopoulos (1975) to 
describe the 2:1 Lindblad resonances in the plane of the 
galaxy but uses the vertical action angle variables. Us- 
ing this model we consider the growth of a boxy/peanut 
shaped bulge via resonance trapping in the vertical inner 
Lindblad resonances caused by the growth of the bar. Fol- 
lowing the growth of the bar, our model predicts that the 
height of the peanut increases as a function of radius until 
it reaches a maximum at the outermost radius which is set 
by the condition | = — 1. Our model naturally accounts 
for the sharp outer edge of observed boxy/peanut shaped 
bulges. The maximum total vertical width of the peanut 
is of order ~ 2\fkv res which we estimate for a typical bar 
could be about two thirds r resi where r res is the radial 
location of the resonance. This model is therefore capa- 
ble of accounting for the observed widths of boxy/peanut 
shaped bulges. 

The scenario explored here suggests that the 
boxy/peanut shaped bulge is grown at the same time as 
the bar forms. We predict that recently formed bars should 
manifest strong boxy/peanut shaped bulges and barred 
galaxies should never be found without boxy/peanut 
shaped bulges. If the boxy/peanut shapes take nearly 
a Hubble time to grow, as found in the N-body simula- 
tions of Combes et al. (1990), then there should be barred 
galaxies which lack boxy/peanut shaped bulges. However 
Bureau & Freeman (1999) considered a control sample of 
edge-on galaxies lacking boxy /peanut bulges. No exam- 
ples of barred galaxies that lack a boxy/peanut shaped 
bulge have yet been identified. Pfcnniger & Friedli (1991) 
found that the boxy/peanut only required 3-4 bar growth 
periods to grow and this could be significantly shorter, 
only a few Gyrs. If this timescalc were more appropri- 
ate then the lack of barred galaxies which do not display 
boxy/peanut shaped bulges is less of a problem. 

One way to probe the formation timescale of the 
boxy/peanut is to search for barred galaxies which have 
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recently formed bars. If boxy/peanuts require many ro- 
tation periods to form then these barred galaxies should 
have normal (non-boxy) bulges. Because of the high rate 
of star formation along the bar Quillen et al. (1997) in- 
ferred that NGC 7582 was a case of a recently formed bar. 
They pointed out that the galaxy bulge did show a strong 
prominent peanut and so argued that boxy/peanut shaped 
bulges are formed quickly after the onset of a bar. Obser- 
vations of edge-on galaxies are therefore consistent with a 
rapid boxy/peanut formation mechanism such as we have 
proposed here. 

Our scenario does not require the galaxy to buckle in 
order to form the boxy /peanut shaped bulge. The fire- 
hose instability scenario for bar-buckling (beautifully il- 
lustrated in 3D N-body simulations by Raha et al. 1991) 
predicts that the galaxy will literally be bent in a U- 
shape while the galaxy buckles, however no examples of 
U-shaped edge-on galaxies have been identified in galaxy 
surveys (with the exception of galaxies which are believed 
to be ram pressure stripped in clusters). It is likely that 
the period of time during which the galaxy is U-shaped is 
short. However galaxy collisions have a similar estimated 
lifetime and are detected in nearby galaxy surveys. We 
note that Friedli & Pfenniger (1990) reported that N-body 
simulations with fixed z-symmetry did grow boxy/peanut 
shaped bulges, though at a slower rate than unconstrained 
simulations. 

We have constructed an analytical framework which can 
be used to predict the height distribution of the peanut as 
a function of radius and the azimuthal angle in the galaxy 
plane. The nature of the dependence is primarily deter- 
mined by the resonance responsible for the vertical exci- 
tations. A 3 dimensional study of a nearby boxy bulge 
(such as that seen in our galaxy) could determine which 
resonance is most likely to be responsible. The shape of 
the boxy/peanut shaped bulge is also sensitive to the de- 
pendence of the gravitational potential on height from the 



mid-plane and on the way the bar grew. Further study 
of boxy/peanut shaped bulges may allow us to place con- 
straints on the vertical shape of the galaxy and on the 
early evolution of the bar. 

The toy model introduced in this paper was designed 
to illustrate the process of resonance capture during the 
growth of a bar. We have argued that the growth of 
the bar itself causes stars to become trapped into verti- 
cal resonances and so lifted out of the plane of the galaxy. 
However, simultaneous growth of a bar and boxy/peanut 
shaped bulge has not yet been observed in three dimen- 
sional N-body simulations of non-bending bars. This could 
in part be because N-body simulations are are typically 
thicker than thin bars such as NGC 7582. Moreover, 
N-body simulations (Pfenniger & Friedli 1991; Friedli & 
Pfenniger 1990; Combes et al. 1990) have shown that the 
boxy/peanut continues to grow well after the bar is sta- 
ble. Since resonance capture is likely to be a natural con- 
sequence of bar growth, it might be useful to reexamine 
N-body simulations and verify if orbits are trapped into 
vertical resonances during the bar growth period. After 
the bar stabilizes, the analytical framework developed here 
might also be applied toward predicting the evolution of 
stellar orbits caused by secular evolution of a bar. For 
example a slow decrease in pattern speed will also trap 
particles in the vertical resonances and cause them to be 
lifted into higher orbits. In this case v decreases with 
time so that phase space moves downwards in Figure 1, 
causing similar phenomenolgy as growing the bar; parti- 
cles trapped into resonances can be lifted higher as the bar 
evolves. 
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APPENDIX 

To put the Hamiltonian into action angle variables to first order we use the generating function 

Fx{z,ff) = -^z 2 tan6»' (23) 
Following this canonical transformation the z part of the Hamiltonian restricted to r = r c 

H z {p z -z) = \{pl + v 2 z 2 ) + ^- (24) 

becomes 

H z {I'-6') = V I> 3 + ¥1 cos 4 (6") (25) 

To put the Hamiltonian into a form that doesn't depend on the angle we must do another canonical transformation. Note 
that cos 4 9 = | cos 46» + ± cos 29 + § . 

We try the following generating function 

F 2 (6',I 3 ) = h- sin26" + 1 sin AO') (26) 

which leads to new variables I3 , 6*3 such that 

/' = 7 3 + A/f Q cos 29' + 1 sin A9'\ (27) 
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Fig. 4. — a) Similar to Fig. 3 except the initial dispersion in jz is a = 0.05. b) Similar to a) except the perturbation was grown in 10 
rotation periods. 
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Fig. 5. — a) For the distribution in action angle variables displayed in Fig. 3 this shows the vertical distribution as a function of azimuthal 
angle for the case m = 2 corresponding to the 1:1 resonance. The x axis shows the azimuthal angle, 9 — Qi,t, in the plane of the galaxy in the 
frame that rotates with the bar and is given in degrees. The y axis shows the z distribution in units of r res . b) Same as a) but for P = 10 
(also shown in Figure 4b). We see from these figures that we expect a maximum peanut width at the radius where j = — 1. 



Fig. 6. — Same as Fig. 5a,b but for the case when 
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To fourth order in J 3 the Hamiltonian becomes 

H z (h;0 s ) = uh + ^ (28) 

which is independent of 6*3 so that I 3 is a conserved quantity and we find that we have successfully transfered the 
Hamiltonian into action angle variables. 

Via a similar canonical transformation the term proportional to (r — r c )z 2 in Equation (5) disappears but the term 
proportional to (r — r c ) 2 z 2 remains yielding the additional term (k 2 — 3f^ 2 ) vI ^ z which we have included in Eqn. (7). 
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